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Abstract 

The physics of nuclear reactions in stellar plasma is reviewed with special 
emphasis on the importance of the velocity distribution of ions. Then the 
properties (density and temperature) of the weak-coupled solar plasma are 
analysed, showing that the ion velocities should deviate from the Maxwellian 
distribution and could be better described by a weakly- nonexstensive (\q — 
1| < 0.02) Tsallis' distribution. We discuss concrete physical frameworks for 
calculating this deviation: the introduction of higher-order corrections to the 
diffusion and friction coefficients in the Fokker-Plank equation, the influence 
of the electric-microfield stochastic distribution on the particle dynamics, a 
velocity correlation function with long-time memory arising from the coupling 
of the collective and individual degrees of freedom. Finally, we study the 
effects of such deviations on stellar nuclear rates, on the solar neutrino fluxes, 
and on the pp neutrino energy spectrum, and analyse the consequences for 
the solar neutrino problem. 



I. INTRODUCTION 

Thermonuclear reactions are of enormous importance for the physics of stars and of our 
Sun. The general formalism that describes these reactions exists since a long time and there 
is a wide consensus about our good understanding of the relevant physics Up-^. 

However, quantitative calculations of specific reaction rates need experimental inputs 
and theoretical assumptions. Cross sections must be known. In the best cases, they can 
be directly measured at the relevant energies; other cross sections can be only measured at 
higher energies and need to be extrapolated; in some situations only theoretical predictions 
exist (for a recent review of the cross sections relevant to the Sun see Ref. M). In addition, 
rates depend on weighted thermal averages of the cross sections; therefore, knowledge about 
the thermal distribution is also needed. 
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Because many nuclear reactions in the stellar burning core proceed by way of quan- 
tum penetration of a high Coulomb barrier, their cross sections grow exponentially with 
energy. Therefore, thermal averages do not probe the average energy of the distribution 
(kT), but its high-energy tail. Consequently, rates are sensitive to a relatively small part of 
the distribution. 

The ion velocity distribution in stellar calculations is always assumed Maxwellian []. In 
fact, it is almost universally accepted as a fact that physical conditions in the solar interior 
(density and temperature) lead to an equilibrium velocity distribution that is Maxwellian. 
In this article we review the assumptions underlying this general statement, discuss why 
they are only approximately true, and argue that the consequent small corrections to the 
distribution have significant effects on the rates. 

The dynamics of the solar plasma is not trivial, since at such densities and temperatures 
there is no clear scale separation between collective and individual degrees of freedom ||. 
The presence of more than just one energy scale {kT) in the relevant range of energies results 
in deviations from the pure exponential behavior, exp (—E/kT), which is determined only 
by kT. We have tried several new approaches to this strongly interacting many-body system 
for which a realistic microscopic calculation does not exist yet. 

In one approach we exploit the knowledge that the distribution in the solar interior 
cannot be too much different from the Maxwellian one and add small corrections (higher- 
order terms in a derivative expansion) to the coefficients of the standard Fokker-Plank 
equation. Tsallis' p|,|T0|| and Druyvenstein-like distributions are immediately generated. 

A second approach || focuses on the electric microfields that have been shown to exists 
in plasmas and tries to link their distribution and the connected effective cross sections with 
deviations from the Maxwellian distribution. Different classes of microfields distributions 
and effective cross sections yield several non-Maxwellian distributions among which we again 
find Tsallis and Druyvenstein-like distributions. 

The third approach has been just started and aims to connect the distribution of collec- 
tive variables |TT| to memory effects and long-time correlations between velocities. There 
should exist solutions compatible with the Tsallis' distribution and/or other non-Maxwellian 
distributions. 

These three approaches are not exhaustive and not necessary alternative. Nevertheless, it 
is suggestive that all of them point in the same direction: the Maxwell-Boltzmann distribution 
of velocity should have small but nonnegligible corrections in the solar plasma and the Tsallis 7 
distribution could provide a better description. 

Nowadays solar modeling seems to have reached a satisfactory stage |I|j3|JT^|. The in- 
clusion in the latest models of higher-order effects, such as the diffusion of heavy elements, 
brings the theoretical predictions in good agreement even with the detailed helioseismological 
data |l3| , p^| . However, there still exist a discrepancy between the solar neutrino experiments 



and the predicted neutrino fluxes |I5|. In this context, there has been a considerable amount 



We consider situations where quantum effects are negligible, e.g., in the solar plasma small 
quantum corrections to the statistics exist for electrons but not for ions. More generally, one could 
also analyse plasma corrections to the standard Fermi or Bose distributions. 
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of work devoted to answering questions such as: how large are the uncertainties of the solar 
model input parameters? Has something been left out of standard solar models? How does 
this affects predictions for the fluxes and the status of solar neutrino problem (SNP) |[16|-F2T 



Therefore, we find extremely important to assess the consequences of the possible small 
deviations from the standard statistics on the solar model and, in particular, on the neutrino 
fluxes. At least two effects should be considerd: changes of the total rates caused by the 
different thermal average and modifications of the shape of the energy spectra. The effect on 
the rates was already considered by Clayton p^j23[| two decades ago. New experimental data, 
better solar models and a considerable better understanding of nonstandard distributions 
and of how they can arise in solar plasmas convinced us of the necessity of reconsider this 
earlier suggestion. Even if experimental detection of small modifications of the neutrino 
spectra might appear still far from being possible, it is useful to have a clear assessment of 
what the signal would be. 

Our paper is organized as follows. In Sec. [TT] we review the physics of subbarrier ther- 



monuclear reactions and how the velocity distribution influences reaction rates. Section [] I 



which is dedicated to the solar plasma, contains the central part of our work: three ap- 
proaches to velocity distribution calculation and the resulting nonstandard distributions. In 
Sec. |TV] we calculate the reaction rates with the modified distributions, while we obtain the 
neutrino spectrum from the pp reaction in the presence of Tsallis' statistics and compare 
it to the standard one in Sec. [V| The solar neutrino problem is analysed in the light of 
our results on nonstandard velocity distributions in Sec. [VI], and Sec. |VH| is reserved to our 
conclusions. 



II. THERMONUCLEAR REACTION RATES 

Reliable calculations of nuclear reaction rates in stellar interiors is fundamental for a 
quantitative understanding of the structure and evolution of stars. In fact, while the overall 
stellar structure is rather robust, changes of some of the rates by few percent can produce 
detectable discrepancies, when precise measurements are possible, e.g., in the case of the 
solar photon and neutrino luminosity, and mechanical eigenfrequencies [El . 



In this Section, we review the basic physics underlying thermonuclear reactions and the 
main ingredients of their calculation M. 



A. Introduction 



Let us consider a gas with n\ particles of type 1 and n 2 particles of type 2 per cubic 
centimeter and relative velocity v ; the reaction rates r (the number of reactions per unit 
volume and unit time) is given by 

r = {I + 5 l2 )~ l n ini {va) , (1) 

where a = a(v) is the nuclear cross section of the reaction. The reaction rate per particle 
pair is defined as 
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poo 

(va) = f(v) av dv , 
Jo 



(2) 



where the particles distribution function f(v) is a local function of the temperature. 

Therefore, the reaction rate per particle pair (va) is determined by the specific cross 
section and by the velocity distribution function of the incoming particles. In general, cross 
sections do not have very strong dependence on the energy, when no energy barrier is present 
and away from resonances. Therefore, most of the contribution to (va) comes from particles 
with energy of the order of kT, and the dependence on the specific form of f(v) is weak. 
We shall see that the situation is very different in the presence of the Coulomb barrier. 



B. Subbarrier reactions 

Most of the nuclear reactions that power the stars are between charged particles. These 
particle must penetrate a Coulomb barrier that is very large in units of kT. The pen- 
etration probability is proportional to the Gamow factor exp [—2tti](E)], where 1](E) = 
ZiZ 2 aJ jjbc 2 1 IE is the Sommerfeld parameter, a is the fine structure constant and \x is the 



reduced mass; this factor can be also written as exp [— y Eq/E], defining the Gamow energy 
Eq = 2fic 2 (ZiZ 2 a7i) 2 . This exponentially small probability makes the cross section grow 
extremely fast with the energy; therefore, one usually defines the astrophysical S factor, 
whose energy dependence is weaker 

a(E) = ffi e -*"»(*> = ^l e ~V^ . ( 3 ) 
E E 

This reaction mechanism has at least two main consequences for the study of stellar 
structure. 

(1) Since the bulk of particles in the stellar plasma have thermal energies of the order of 
kT, which is far below the Coulomb barrier, only a small number of particles in the high- 
energy tail of the distribution has a chance of reacting: this high-energy tail plays a crucial 
role for the reaction rates, which becomes very sensitive to small changes of the distribution. 

(2) The cross sections of some of the reactions, e.g., weak reactions with Coulomb barrier 
such as p + p — > d + e + + v, are extremely small at the low energies relevant to stellar physics 
(of the order of tens of keV for reactions in the Sun). Therefore, some of them have only 
been calculated theoretically, others have only been measured at higher energies and then 
extrapolated to the lower thermal energies of the stellar interiors. 

Both these two facts lead to considerable uncertainties in the reaction rates. The second 
point has been extensively discussed, and it is common opinion that the present standard 



theory of stars, and in particular of the Sun pTlfl] , has already taken properly into account 



the uncertainties in the cross sections. On the contrary, the first point, the possibility that 
small changes of the particle energy distribution could strongly affect reaction rates, has 
received little attention [p2j25,p6,5|, or dismissed on the ground that the energy distribution 



in stars is allegedly well-know 0]. In this paper, we shall discuss mainly this point, the 
energy distribution in stellar plasmas and its implication for nuclear rates. 



4 



C. Maxwell distribution and the Gamow peak 



Normal stellar matter, such as the one in the Sun, is nondegenerate, i.e., quantum effects 
are small (in fact, they are small for electrons and completely negligible for ions), it is 
nonrelativistic, and it is in good thermodynamical equilibrium. On this ground, the particle 
velocity distribution is almost universally taken to be a Maxwell-Boltzmann distribution, 
without much questioning. 

However, derivations of the ubiquitous Maxwell-Boltzmann distribution are based on 
several assumptions ||. In a kinetical approach, one assumes (1) that the collision time 
be much smaller than the mean time between collisions, (2) that the interaction be suffi- 
ciently local, (3) that the velocities of two particles at the same point are not correlated 
(Boltzmann's Stosszahlansatz) , and (4) that energy is locally conserved when using only the 
degrees of freedom of the colliding particles (no significant amount of energy is transferred 
to collective variables and fields). In the equilibrium statistical mechanics approach, one 
uses the assumption that the velocity probabilities of different particles are independent, 
corresponding to (3), and that the total energy of the system could be expressed as the sum 
of a term quadratic in the momentum of the particle and independent of the other variables, 
and a term independent of momentum, but if (1) and (2) are not valid the resulting effective 
two-body interaction is not local and depends on the momentum and energy of the particles. 
Finally, even when the one-particle distribution is Maxwellian, additional assumptions about 
correlations between particles are necessary to deduce that the relative- velocity distribution, 
which is the relevant quantity for rate calculations, is also Maxwellian. 

At least in one limit the MB distribution can be rigorously derived: systems that are 
dilute in the appropriate variables, whose residual interaction is small compared to the one- 
body energies. In spite of the fact that the effects of the residual interaction cannot be 
neglected, as a good first approximation the solar interior can be studied in this dilute limit; 
therefore, it is reasonable to suppose that the velocity distribution in the Sun is not too far 
from the Maxwellian one. This fact is consistent with the many successes of the standard 
treatment and suggests that we start discussing the standard results and improve them later. 

In the ordinary treatment, the single particle energy distribution for protons and other 
ions is taken as 



w*) = £^«-'». (4) 

When this distribution and a cross section of the form of Eq. (|3]) are inserted in Eq. (0) , the 
resulting integrand goes to zero both at large energies, because of the exponentially small 
number of particles, see Eq. (|J), and also at small energies, because of the exponentially 
small probability of barrier penetration, see Eq. @. In fact, the integrand has a maximum 
at the temperature-dependent energy 

^{Mf.f^ TkT , (5) 

which it is called the most effective energy [pj, since most of the particles that react have en- 
ergies close to E ] we have defined the adimensional parameter r = E /kT = [E G / (4/cT)] 1 / 3 . 
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Figure [l| gives a pictorial demonstration of how the Gamow peak originates. The ex- 
ponentially decreasing function (solid curve labeled "Maxwell") is the Maxwellian factor 
exp (—E/kT) multiplied times 400 to emphasize the tail of the distribution. The rapidly 
growing function (dotted curve labeled "Penetration factor") is the factor exp — <J Eg / E 
with the choice of Eg = 4000/cT, which corresponds to a most effective energy Eq = 
(EG(kT) 2 /4) 1 / 3 = lOfcT; it has been arbitrarily normalized such that it is equal to 400 
at E — 20kT. The product of the two functions yields the Gamow peak (solid curve labeled 
Gamow), which has been normalized to one at its maximum (E — Eq — 10kT). It is impor- 
tant to notice that the area under the Maxwellian curve for energies larger than about QkT 
(where one finds most of the contribution to the peak) is less then 0.3% of the total area. 

In this framework, one performs the integral in Eq. (0) using the saddle-point asymptotic 
expansion around the maximum Eq of the Gamow peak and finds that the reaction rate per 
particle pair is, apart for small calculable corrections, 

W12 = n 8 — S(E )T* 2 e-^ . ( 6 ) 
V3/ii2C7raZiZ 2 

Table | reports the values of E in units of kT (r) for several reactions. We notice that 
Eq can be much larger than kT and, therefore, only a very small number of particles in the 
far tail of the distribution contributes to the rate, e.g., for the 3 He + 3 He reaction the most 
effective energy is about YJkT: there are less than about 40 particles out of a million that 
have energies so large or larger. It is not sufficient anymore to know that the Maxwellian 
distribution is good approximation, we must be sure that there are no corrections to a very 
high accuracy. 



III. VELOCITY DISTRIBUTIONS IN NON-IDEAL PLASMAS 



As discussed in Sec. |Q[ several subtle points must be assumed in the derivation of the 
often-taken-for-granted Maxwellian velocity distribution. Likewise, one often assumes that 
the solar core could be treated as an ideal (Debye) plasma. However, there are physical con- 
ditions and/or specific applications that needs higher accuracy for which becomes necessary 
to take into account modifications of the standard plasma theory. 

In this Section we discuss how the physics of non-ideal plasmas, and in particular of 
the solar interior, can result in equilibrium velocity distributions that deviate from the 
"standard" one. In addition, we present specific physical frameworks where non-Maxwellian 
distributions arise and microscopic mechanisms that allow a reliable estimation of the size 
of deviations. 



A. Ideal and non-ideal plasma 

In literature, a plasma is characterized by the value of the plasma parameter Y 

a kl 
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where a = n" 1 / 3 is of the order of the interparticle average distance (n is the average density). 
The plasma parameter is a measure of the ratio of the mean (Coulomb) potential energy 
and the mean kinetic (thermal) energy. 

Depending on the value of the plasma parameter, we can distinguish three regimes 
that are characterized by different effective interactions and require different theoretical 
approaches. 

• T < 1. The plasma is described by the Debye-Huckel mean-field theory as a dilute 
weakly interacting gas. The screening Debye length 



/ kT 

is much greater than the average interparticle distance a, hence there is a large number 
of particles in the Debye sphere (N D = (A-k/3)R 3 d ). Collective degrees of freedom are 
present (plasma waves), but they are weakly coupled to the individual degrees of 
freedom (ions and electrons) and, therefore, do not affect their distribution. Binary 
collisions through screened forces produce the standard velocity distribution. 

T pa 0.1. The mean Coulomb energy potential is not much smaller of the thermal 
kinetic energy and the screening length Rp ~ a. It is not possible to clearly separate 
individual and collective degrees of freedom. The presence of at least two different 
scales of energies of the same rough size produces deviations from the standard statis- 
tics which describe the system in terms of a single scale, kT. 

r > 1. This is a high-density/low-temperature plasma where the Coulomb interaction 
and quantum effects start to dominate and determine the structure of the system. 



B. The solar interior 

In the solar interior the plasma parameter is T pa 0.1; therefore, the solar core is a weakly 
non-ideal plasma where the Debye-Huckel conditions are only approximately verified. Similar 
behaviors are expected in other astrophysical systems characterized by plasma parameters in 
the range 0.1 < T < 1; examples are brown dwarfs, the Jupiter core and stellar atmospheres. 

Studies of systems with this intermediate values of T are the most difficult; it is possible 
to use/combine several different approaches none of which, however, is completely justified. 

The reaction time necessary to build up screening after a hard collision can be estimated 
from the inverse solar plasma frequency t pi = uj~^ = \Jm/ (Anne 2 ) pa 10~ 17 sec, and it is 
comparable to the collision time t co u = (avn)^ 1 pa 10~ 17 sec. Therefore, several collisions 
are likely necessary before the particle looses memory of the initial state and the scattering 
process can not be considered Markovian. In addition, screening starts to become dynamical: 
the time necessary to build up again the screening after hard collisions is not negligible any 
more. 
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C. Three roads towards nonstandard distributions 



For concreteness we consider three possible approaches to weakly-noninteracting plas- 
mas that yield nonstandard distributions. These approaches are not alternative, but they 
could perhaps be seen as different and partial descriptions of the same complicated physical 
problem. The main purpose of their presentation is to give concrete examples of how non- 
standard statistics arises. However, they are not by any mean the last word on the topic: we 
must still develop their full potentiality and alternative approaches must also be pursued. 



1. Fokker-Plank 

In the Fokker-Plank context it is possible to introduce corrections (i) to the lowest order 
friction coefficient J(v), and (ii) to the lowest order diffusion coefficient D(v) |||7]>0]- Similar 
corrections have already been shown to exist in hydrodynamic systems. 

We assume that the system is not too far from the standard regime that leads to the 
MB distribution, so that an expansion starting from the usual formalism makes sense. The 
Fokker-Planck equation, given in the Landau form, is 

|/(*. «) = I; ('MM, v) + ^D(v)f(t, t,)) , (9) 

where fit, v) is the distribution probability of particles with velocity v at time t and J(v) and 
D(v) are the dynamical friction and diffusion coefficients. The stationary distributions are 
the asymptotic solutions of the above equation. To lowest order J{y) — v/r and D{v) = e/r, 
where the constant r > has dimension of time (m/r is the friction constant) and \fe has 
dimension of a velocity (e = kT/m for Brownian motion). At equilibrium one obtains the 
well-known Maxwellian distribution 

lim /(M )~exp{-g=exp{-^} . (10) 

We can generalize the standard Brownian kinetics considering the expressions of the 
quantities J(v) and D[v) to the next order in the velocity variable: J{v) — v/r (1 + Piv 2 /e) 
and D(v) = e/r (1 + 7if 2 /e); these higher derivative terms can be interpreted as signals of 
nonlocality in the Fokker-Planck equation. 

If Pi = and 71 7^ we find the Tsallis' distribution 



/(*) 



1 + 07-1) 



mv 2 



2kT 



e 



1 + (q - 1 



mv 2 
2kT 



(11) 



where q — 1 = 271/(271 + 1), is the Heaviside step-function, and kT/m = e(2 — q). When 
the characteristic parameter q is smaller than 1 (—1/2 < 71 < 0), this distribution has a 
upper cut-off: mv 2 /2 < kT / (1 — q) (the tail is depleted). The distribution correctly reduces 
to the exponential Maxwell-Boltzmann distribution in the limit q — > 1 (71 — > 0). When 
the parameter q is greater than 1 (71 > 0), there is no cut-off and the (power-law) decay is 
slower than exponential (the tail is enhanced). 
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If Pi 7^ and 71 = 0, we find a Druyvenstein-like distribution: 

/M~exp{4-ftg) 2 }, (12) 

which has also the functional form suggested by Clayton to parameterize a small devi- 
ation (depletion) from the Maxwellian statistics. 



2. Random fields 



Each particle is affected by the total electric field distribution due to the other charges 
in the plasma. If external fields and large scale internal fields due to collective modes 
are subtracted from the total electric field, the single particle see the remaining field as 
a relative small random component. The density of these random electric microfields has 
been studied [27J28| and it is often expressed in terms of the dimensionless parameter F as 
(S 2 ) = (Fe/a 2 ) 2 , where 4ira 3 /3 — 1/n; the distribution of F in plasmas depends on the 
value of T ||28|| . These microfields have in general long time correlations, and can generate 
anomalous diffusion. 

The total (micro)field can be decomposed into three main components. 

(i) A slow-varying (relative to the collision time) component due to collective plasma os- 
cillations, which the particle sees as an almost constant external mean field £ over several 
collisions. 

(ii) A fast random component due to particles within a few Debye radii, whose effect can 
be described by an elastic diffusive cross section <j\ ~ v . When only this cross section is 
present, the distribution remains Maxwellian even in presence of the slow mean field 8. 

(iii) A short-range two-body strong Coulomb effective interaction, that can be described by 
the ion sphere model [29|. The strict enforcement of this model as implemented by Ichi- 
maru [^9| yields the elastic cross section cr = 2na 2 a 2 , where a is the interparticle distance, 
a an adimensional parameter whose order of magnitude can be inferred from the parameter 
F that characterizes the microfields, F ~ or 2 . Since F 2 ~ 3/T « 40, for T = 0.07, one can 
estimate 0.4 < a < 1. In the present model, it is this component of the electric field that 
turns out to be mainly responsible of the correction factor exp [— 5(E/kT) 2 ]. 

In this framework, the stationary solution of the kinetic equation valid for small devia- 
tions from the MB distribution can be shown to be the Druyvenstein distribution: 



f(E) ~ exp 



-<P 



E 



'■JO' 



(13) 



where 



1 + V 



-H 



(14) 



'3<o?> 




3<o?) 



(15) 
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k = 2m a mb/(m a + mb) 2 is the elastic energy-transfer coefficient between two particles a and 
b, and S is proportional to the square of the ratio of the energy densities of the electric field 
and of the thermal motion. 

In the small correction limit, relevant to the solar interior, = 1 and the 5 parameter is: 

|5|^-f^ = 12« 4 r 2 «l. (16) 

From Eq. (0), we see that the presence of electric microfields implies a deviation from 
the Maxwell-Boltzmann distribution and that the entity of this deviation depends on the 
value of the plasma parameter. As already mentioned, in the solar core T m 0.1, hence, 
the order of magnitude of the deviation parameter 5 is about 0.01 and in this intermediate 
region long-range and memory effects take place. 



3. Memory effects and collective variables 

As already discussed, in the solar core collective effects have time scales comparable to 
the average time between collisions (e.g., compare the inverse plasma frequency with the 
average collision time lu^ 1 ~ t co u) and it is not possible a description that separates the 
collective and the individual degrees of freedom. 

A rigorous microscopic approach would imply the resolution of the dynamical equations 
of motion with a Hamiltonian that explicitly contains collective and individual degrees of 
freedom and their mutual interactions. Our scope is more modest: we want to give a concrete 
example where memory effects (long-time correlations) are important. 



The authors of Ref. 30, IT! introduce the collective variables 



N 

Pk(t) = ^2esxp{-ikXj(t)) , (17) 

3 

and describe these collective variables as harmonic oscillator variables weakly coupled to 
the individual degrees of freedom, which, therefore, act as a thermostat for the collective 
variables. 

Generalizing this approach, we write down a generalized Langevin equation for the col- 
lective variables 

Pk = -ujl Pk + [ K(t-r)p k (r)dT + F{t) , (18) 
Jo 

where K(t — r) is the memory friction kernel, which takes into account the long-time tail 
of the correlations in the solar interior, and F(t) is the stochastic fluctuating thermal force 
due to the interactions with the individual degrees of freedom. In weakly non-ideal plasmas, 
the longitudinal waves have the following dispersion law 

u 2 k = u 2 pl (i + m 2 D k 2 ) . (19) 

The two functions K(t — r) and F(t) are not independent but related by the second 
fluctuation-dissipation theorem 
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2kT 

(FfaWh)) = K{t x -t 2 ) . (20) 

In the limit of negligible memory effects, the memory kernel becomes a 5-function, that is, 
K(t) = 7#(t), where 7 is the effective friction coefficient, which can be written in terms of 
the Landau (7^) and collisional (j c ) damping JTTJ . 



The explicit form of the memory kernel depends on the specific system and, in principle, 
should be deduced from a microscopic calculation. At the moment, we have been able to 
verify that, in the framework of Tsallis nonextensive statistics, the memory kernel can be 
written, for q sa 1, as 

K(t)= 1 6(t) + (l-q) 1 2 ^l + 2 1 t+tPj . (21) 

The above expression is consistent with the prescription of the nonextensive statistics 
and implies anomalous subdiffusion for q < 1 and superdiffusion for q > 1. A more complete 
analysis of this possibility is still under way; it is also important to study the effects of such 
a memory kernel for the collective modes on the individual degrees of freedom. 

At the moment, we can speculate that such correlations among the collective modes 
could lead to a long-time asymptotical behavior of the velocity-correlation of the ions of the 
kind (v(0)v(t)) ~ t~( 1+1 \ If 7 > 1, i.e., the correlation decays sufficiently fast, the diffusion 
process is no qualitative different from the delta- function case: (x 2 (t)) ~ t. In this case, 
one can show PHJ52]] that if < 7 < 1 and < / (v(0)v(t)) < 00, the standard distribution 
of velocity is still valid. However, if < 7 < 1 and J (v(0)v(t)) = (or very small), or 
if — 1 < 7 < 0, the diffusion is anomalous (x 2 (t)) ~ t 1+1 (~ tlogt, if 7 = 0). Indeed, 
Tsallis [|33] shows that the generalized entropy quite naturally generate both anomalous 
diffusion ((x 2 (t)) ~ t 1+7 ) and the non-Maxwellian probability distribution for the velocities 
ofEq. (P). 



D. Summary 

In this Section we have presented theoretical arguments that the velocity distribution in 
the solar core should deviate from the standard one and, in particular, that it could follow 
Tsallis' distribution p|JT0|] shown in Eq. fllTl) . We have also estimated the deviation from the 



standard statistics from the known microfield distribution. If we use the parameterization 



suggested by Clayton 0,|3|] (Druyvenstein distribution) 

f{E) ~ (kT)~ 3 / 2 e -E/kT-S{E/kTf ? (22) 

S should be of the order of 0.01. Since in the limit q — > 1 the Tsallis distribution can be 
asymptotically described as a Druyvenstein with 5 = (1 — q)/2, we have also estimated that 
q should be a few percent different from 1. 

In the next section, we shall demonstrate how even such small deviations from the 
Maxwell-Boltzmann distribution can be very important for solar physics. 
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IV. REACTION RATES AND MODIFIED THERMAL DISTRIBUTION 



Both from general considerations about the successes of the standard approach and from 
the estimates and calculations that we have shown, it should be clear that the deviations 
from the Maxwellian distribution in the Sun are small. Therefore, it is completely general 
to make an asymptotic expansion of the following kind: 

f{E) ~ (kT)- 3 / 2 e -e/kT-&(E/kTf ; (23) 

where we have disregarded terms with powers higher than (E/kT) 2 in the exponent and 
all other power corrections outside the exponential, apart the ones (not shown) needed to 
correctly normalize the distribution. This result becomes the more accurate the more 5 is 
small. Note that the exponential can not be expanded. 



This same kind of parameterization was already considered by Clayton [23]. One can 
easily convince oneself that the distributions previously discussed can be put in this form 
in the limit of small deviations. We have seen that the the Tsallis' distribution plays a 
special role: it can be also approximated to first order in (1 — q) by Clayton's form with 
5 = (1 — q)/2 and a renormalized temperature T' = T + T(l — q). In Sec. II we performed 
a saddle point expansion of the integral over the velocity distribution with the introduction 
of the most effective energy, E , and approximating the integrand around the maximum E 
with a Gaussian function (Gamow peak). The same kind of asymptotic expansion can be 
repeated for this modified distribution yielding the following analytical expression for the 
change of the rate: 

^ = e -^o } /^ 2 ee e~** , (24) 
where E is the most effective energy of Eq. (H) 



E MM T c ^/ 3 



, 5 ,64 [Z\Zl , 1 \ ^ , (25) 
kT V A 1 + A 2 T J y ' 

which depends on the reaction i, through the charges Z and weights A of the ions, and on 
the relevant average temperature T (T c = 1.36 keV is the temperature at the center of the 
Sun), and where (v<Ji)o is the expression in Eq. (|j). 

We shall use the expression in Eq. ( |2~3] ) to discuss the effects of modified statistics on 
the solar thermonuclear rates: it has a sufficient numerical accuracy for this situation and it 
is physically more transparent. However, we perform the integral over the distribution nu- 
merically in physical situations where stronger deviations from the Maxwellian distribution 
are expected, and when we want to check whether the asymptotic expansion is sufficiently 
accurate. 

The dramatic effect of small deviations from the Maxwellian distribution on the rate can 
be appreciated by looking at Fig. |I[ The steep decreasing dashed curve shows the high- 
energy tail of the Tsallis' distribution [1 + (q — 1)E /kT] 1 ^ 1 "^ with q = 0.98 compared with 
the corresponding Maxwellian distribution (solid curve). The tails of the two distributions 
shown in Fig. [I] contain about 0.30% (Maxwell) and 0.26% (Tsallis) of their total area and 
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have been multiplied times a factor 400 to make their small difference more visible. The 
product of each of these two very similar distributions times the same penetration factor 
(dotted curve) yields two very different Gamow peaks. The peak corresponding to the Tsallis 
distribution (dashed curve) is not only shifted at lower energy by more than one unit of kT 
compared to the Maxwellian peak (solid curve), but it is also much lower: a tiny change of 
the distribution tail (A x 10~ 4 less probability in the tail) makes the probability in the peak 
smaller by a factor of about 0. 31 

V. THERMAL EFFECTS ON THE NEUTRINO ENERGY SPECTRUM 

Because of their small interaction cross section, neutrinos that reach the earth have not 
interacted with the solar matter and carry direct information from the solar core. One 
might hope that their energy spectrum could tell us something about the particle velocity 
distribution in the Sun. 

The neutrino energy spectrum for a given reaction depends on the total energy available 
in the center of mass to the final products, total energy that is the sum of the thermal energy 
of the incoming particles and of the energy released by the reaction, and on the velocity of 
the center of mass that produces Doppler broadening. As a specific example, we shall use 
the pp fusion reaction, since its neutrino spectrum is easily calculable from phase space and 
since it releases a relatively low excess energy leaving more space for thermal effects. 

The reaction 

p + p-+d + e + + v e (26) 

releases an excess energy of Q = 2m p -m d ~m e = 2x938.272-1875.613-0.511 = 0.420 MeV, 
which is shared among the kinetic energies of the outgoing particles. Note the the total 
energy released by the reaction is 1.442 MeV, since the positron annihilation gives additional 
1.022 MeV; this energy, however, does not contribute to the kinetic energy of the neutrino 
(this should be contrasted to the more rare reaction p + p + e — > d + v e ). 

Since there are three particles in the final state, the neutrino energy in the center of 
mass has a known probability distribution. For simplicity, we report only the result that 
disregards the recoil energy (a relative correction of the order of 2 x 10~ 4 ) 

P K {E V ) ~E 2 u (K + Q + m e - E u )^j{K + Q - E V ){K + Q + 2m e - E~) x F(E V ) 

= El E e p e x F(E U ) . (27) 

The neutrino energy goes from zero to to the maximal energy 

K - = l K + Q ) {l- j ^ jR )«K + Q, (28) 

where K is the relative kinetic energy of the incoming protons, the approximate equality 
disregards the recoil energy consistently with Eq. (^7|), and where F is the Fermi function 
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F(E U ) = Ee/Pe 

exp [2-ira E e /p e ] — 1 



exp 



2na(K + Q + m e -E l/ )[(K + Q- E V ){K + Q + 2m e - E,)}^ 2 
=T72] 7 ' l 29 J 



2na(K + Q + m e -E u )[(K + Q- E U )(K + Q + 2m e - E v )\ 



1 



which takes into account the Coulomb repulsion of positron and deuteron in the final state. 

The neutrino energy spectrum is obtained by convoluting this result with the distribution 
for K: S(E U ) ~ / dK Pk{E v ) f {K) . In Fig. || we show the energy distribution of the 
neutrinos from the pp reaction for three cases: (1) T = (dotted curve), (2) T = 1.36 keV 
with a Maxwellian distribution (solid curve) and (3) T = 1.36 keV with a q — 1.2 Tsallis' 
distribution (dashed curve). 

Note that relative large value of \q — 1| ~ 0.1 are necessary in order that the effect on 
the spectrum of having a different distribution be of the same order of magnitude of the 
standard thermal effect. 

In principle, precise measurements of the neutrino spectrum give information on the 
thermal distribution of the proton kinetic energies f{K). Let us estimate the necessary 
accuracy. The central temperature of the Sun is about kT c = 1.36 keV and the most 
effective energy Eq (maximum of the Gamow peak) is about five times as big for the pp 
reaction (see r in Table H); therefore, thermal effects give a relative change of the maximal 
energy about 5 ■ 1.36/420 ~ 0.02 and, in general, modify the shape of the spectrum by 
a few percent. The center of mass motion gives also contributions of order kT. Spectral 
measurements with accuracy better than 1% start to be sensitive to the particle thermal 
motion and can measure the solar internal temperature. However, measurements that are 
two orders of magnitude more accurate are necessary to measure the temperature with a few 
percent error and/or discriminate between the Maxwellian distribution and those slightly 
distorted distributions (\q — 1| ~ 0.01 or less) that we shall consider in the next paragraph. 
Such kind of experiments are not expected to be feasible in the near future. 



VI. SOLAR NEUTRINO EXPERIMENTS AND THE SOLAR NEUTRINO 

PROBLEM 

A. The solar neutrino problem 

The solar neutrino problem is one of the most interesting long standing puzzle of the 
modern physics. The combined results from the solar neutrino experiments (Homestake, 
GALLEX, SAGE, Kamiokande and SuperKamiokande) cannot be reconciled with the predic- 
tions of the standard solar models (SSM). Since SSMs have been very successful in predicting 
the stellar structure and have given excellent descriptions of measurements as detailed and 
accurate as the helioseismic ones, the neutrino experiments have suggested that the mini- 
mal standard electroweak model should be extended and include small neutrino masses and 
lepton flavor nonconservation |T5|j2l[] . Neutrino oscillation theory would have far-reaching 
consequences for both particle physics and cosmology: therefore, it is of great importance to 
ask whether the solar neutrino problem can be solved, or at least alleviated, in the framework 
of the conventional physics . 
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The problem can be more precisely appreciated with the help of Fig. [3] and Table pi . 
Figure ^| show the energy spectrum of the main neutrino fluxes coming out of the Sun. The 
same figure shows also the part of the spectrum that contributes to each of the three present 
neutrino measures. The energy dependence of each single flux depends practically only on 
nuclear physics and it is independent of the SSM (see the discussion about the pp spectrum 
in the Sec. |V|). However, the integrated fluxes, which are reported in the first column of 
Table |TJ, and, therefore, the relative weight of the different components of the spectrum, 
depend on the SSM. 

One should compare these predictions with the experimental data (last column of Ta- 
ble imp. 

Only one experiment, the one performed by Kamiokande and SuperKamiokande, (de- 
tection of the Cherencov light emitted in water by electrons that are elastically scattered 
by solar neutrinos, for the latest updated results see Ref. fl35| ) measures neutrinos from a 
single solar reaction ( 8 B — > 2a + e + + u e , the so-called boron neutrinos), since its thresh- 
old (see Fig. |3|) is too high to see any of the other lower energy neutrinos. Therefore, its 
determination of the boron flux is solar- model independent; its present result corresponds 
to less than one half of the predicted flux. This experiment has measured also the energy 
spectrum of the boron neutrinos confirming the expected shape apart for same discrepancy 
for energies greater than about 14 MeV, discrepancy that could be interpreted as evidence 



for oscillations [36 |. 



The other two kinds of (radiochemical) experiments [p?| -|39[| have lower thresholds and, 



since they do not have energy information, measure a combination of several neutrino fluxes. 
Their interpretation in terms of individual fluxes partially depends on the relative weight of 
the fluxes in the solar model. However, they also show lower signals than predicted. For a 
review on the SNP see Refs. [pip! . 



B. Standard attempts of astrophysical solutions 

The SSM flux predictions have ranges of variability, whose magnitudes mostly depend on 
their importance in the energy production mechanism. The pp flux is almost independent of 
the SSM (a few percent uncertainty), the 7 Be is fairly stable (about 10% uncertainty), while 
the boron flux, which is produced in a marginal chain of the energy production, has a larger 
(more than 30%) uncertainty. The ranges of predictions of the SSM fluxes come mainly from 
the experimental/theoretical uncertainties in the input parameters (nuclear cross sections, 
photon opacities and initial contents). 

The have been numerous attempts of reconciling solar models and neutrinos experiments 
by allowing all theoretical inputs to vary within and also considerably outside their uncer- 
tainties. Others have assumed arbitrary changes of the internal structure of the models 
(scaling of the temperature profile, fast element mixing, etc.). This approach to the SNP is 
called the "astrophysical solution" (references could be found in Ref. pip. 



The most radical astrophysical solution consists in leaving the neutrino fluxes as free 
parameters (no solar model) with the only constraint that the produced energy should 
match the solar luminosity. Nevertheless, there remains a discrepancy at about the 3- 
4 a level (the central predictions of the SSM are more than lOc's from the experimental 
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results) IDMOUOMIia. 



C. Non standard velocity distributions 

A different velocity distribution yields a different (va), as discussed in Sec. [TV], and, 
therefore, a different solar model and different neutrino fluxes. 

In principle, the same variation of (va) that is obtained by a different velocity distri- 
bution can also be obtained by an appropriate change of the cross section, for instance by 
means of the astrophysical factor S in Eq. (H). Therefore, allowing a modified velocity dis- 
tribution does not span a different set of neutrino fluxes than the arbitrary variation of all 
the astrophysical factors. In particular, one cannot hope to achieve a better agreement with 
the data than the one achieved by leaving the fluxes themselves as free parameters. 

In practice, however, arbitrary large variations of the astrophysical factors contrast with 
our theoretical understanding of the nuclear reactions and often with experimental mea- 
surements. Allowing generous, but not arbitrary, ranges for the input parameters, implies a 
much more restricted values of the possible neutrino fluxes. 

On the contrary, small changes of the velocity distribution can give huge changes of 
(va). Therefore, changing the velocity distribution has several motivations. (1) When 
one is able to calculated the velocity distribution from the microscopic physics, it gives a 
predictions for (va). (2) Even when one can only estimate the size of the deviation, it gives 
a physical mechanism to justify large and correlated changes of (va). (3) In any case, a 
given uncertainty in form of the distribution yields calculable greater ranges of solar model 
predictions compared to the ones of the SSM, which assumes the distribution function known 
to high accuracy. 

Changing (va) for the zth reaction will affect the whole solar model and, in general, all 
fluxes will change. We use Eq. ([24]) as a general parameterization of the rate change for 
small deviations from the standard distribution. The effects of modifying the rates on the 
fluxes can be estimated by using the power-law dependences 

_1 A(^>oJ ' (30) 

for the fluxes j = 7 Be, 8 B, 13 N and 15 0, and using the solar luminosity constraint [21J to 
determine the pp flux, R pp = 1 + 0.087 x (1 - R Be ) +0.010 x (1 - R N ) +0.009 x (1 - R Q ), and 
keeping fixed the ratio £ = Q pep /&p P = 2.36 x 10~ 3 . The exponents = <91n$j/<91n(t><7j) 
[see Table UD have been taken from Ref. 21], where it is also discussed why solar models 



depend only on the combination (va)u/ J (va)^ and why it is a good approximation to keep 
the ratio £ constant. 

A direct microscopic calculation would determine 5 that could be different for every 
reaction (5 — > Si). The energy distribution can be influenced by the specific properties of 
the ion (charge and mass) and by the different conditions of the environment in those parts of 
the Sun where each of the reactions mostly takes place. However, such a direct calculation is 
not simple and it does not exist for the solar interior. Therefore, for the purpose of estimating 
the effect of nonstandard distributions, we consider two simple models and consider <5('s) as 
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free parameter (s). The first model assumes the same deviation 8 for all distributions, the 
second model assumes that only the p + 7 Be and 3 He + 4 He relative energy distributions 
are nonstandard and introduces Sr^ and 8^ to parameterize their deviations. 
In the first case, one finds by substituting Eq. fl24]) into Eq. (|30"D that 



(31) 



where (3j = J2i ®-jC1% are reported in Table |TJ. This dependence of the fluxes on 8 is in 



good agreement with Clayton's numerical calculation [p3fl . Using the model of Ref. [12] as 
reference model and the experimental results up to the end of 1997 (see Table we obtain 
the best fit for 8 = 0.005 with a \ 2 = 35. 

In the second case, we proceed similarly, but we use 8^ for the reaction p + 7 Be 
and 5(34) for the reaction 3 He + 4 He: the corresponding [3^ = J2i=i7 an d fif 4 ^ = 
^=34^,34(734 — 733) are also reported in Table |TJ. As shown in Table fTJ the best fit is 
obtained for 8^ = —0.018 (negative 8 corresponds to an enhanced tail, q > 1 in Tsallis' 
distribution) and 5(34) = 0.030 with a x 2 = 20. 

As expected from the discussion in the first parts of this Section, this result is a solution 
to the SNP, in the sense of providing a model to fit the experimental results within one (a 
few) sigma. However, we have shown that deviations from standard statistics corresponding 
to values of 8 of about 1% can change the neutrino fluxes of factors comparable to those that 
constitute the SNP. Perhaps, the actual values of the neutrino fluxes coming out of the Sun 
could result from the interplay of several mechanisms that are disregarded in the standard 
picture |34 . 



Moreover, it is clear that the uncertainties of the neutrino fluxes are considerably under- 
estimated by not considering the possibility of non-extensive distributions. 



D. Helioseismic constraints 

Helioseismology provides very detailed and precise information on the solar structure. 
The extremely precise measurements of a tremendous number of frequencies give the possi- 
bility of extracting values of the sound speed even near the solar core where the energy is 
generated. In addition, several properties of the convective envelope are accurately deter- 
mined. 

It is important to verify that nonstandard distributions do not contrast with such data. 
This study has been done for the velocity distribution of proton jEJ: if 8 is within with 
following limits 

- 4.9 x 10~ 3 < 8 PP < 2.3 x 10~ 3 (32) 

there is no incompatibility between helioseismic data and nonstandard statistics. 

As we have seen, even such small values of 8 have very important consequences for the 
neutrino fluxes; in addition, this limit does not automatically apply to the distribution of 
other ions. 
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VII. CONCLUSION 



The solar core, where its energy is produced, is a weakly nonideal plasma. 

Many approaches are possible to such systems. The standard Debye-Hiickel theory is 
very successful, but it is not sufficient to give an accurate description when one needs to 
calculate the energy (velocity) distribution function to a high accuracy. 

We have tried several new approaches to this complex problem. In particular, we have 
considered corrections to the Fokker-Plank equations, the known stochastic distribution 
of electric microfields, and memory effects arising from the interaction of individual and 
collective variables. 

All our attempts indicate that the velocity distribution of ions in the plasma could be 
different from the Maxwellian one, and should be well described by a Tsallis' distribution 
slightly nonextensive (\q — 1| ~ 0.01). 

Such small deviations from the standard statistics produce effects on the energy depen- 
dence of the neutrino spectra or on the helioseismic observables that are not in contrast with 
present data. 

However, even such small deviations of the energy distribution produce dramatic effects 
on those nuclear rates whose main contributions come from the high-energy tail of the 
distribution, as it is best exemplified by Fig. [I]. 

This theoretical possibility enlarges the range of predictions for the solar neutrino fluxes 
and, while it is not sufficient to solve the solar neutrino problem, can make it somewhat less 
dramatic. 
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TABLES 



TABLE I. Most effective energies for thermonuclear reactions and exponents 7 that character- 
ize the change of the thermal average (va) to the leading order in 5, when the energy distribution 
changes by a factor exp{— 5(E/kT) 2 }: (va)s = (va)o exp{— £7}. 

reaction r = Eo/kT 7 = t 2 



(va)n: p + p 2 H + e+ + v 4.8 23 

(va) 17 : p + 7 Be^ 8 B + 7 13.8 190 

(va) 33 : 3 He + 3 He^a + 2p 16.8 281 

{va) u : 3 He + 4 He -> 7 Be + 7 17.4 303 

(va)i M : p + 14 N^ 15 + 7 20.2 407 



TABLE II. The first four rows show = d]n$>j/d\n{vcr)i, the logarithmic partial derivative 
of neutrino fluxes with respect to the parameter shown at the left of the row. These numbers 



are discussed in Ref. [21|. The last three rows show (3j, (3^ e and /^j 1 * 5 , the logarithmic partial 



derivative of the fluxes with respect to the parameters <J's; as discussed in the text, they are linear 
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TABLE III. The first three columns show the predicted fluxes, and the predicted gallium and 



chlorine signals in the SSM [12] and in the two models with nonstandard distribution described 
in the text. The last column shows the present experimental results. For the three models is also 
given the \ 2 resulting by the comparison with the experimental data. 









Models 
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2.55 ±0.21 a 
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137.0 


100 




97 


75 ±5 b 


chlorine 


9.3 


2.84 




3.34 


2.54 ±0.20 c 


x 2 


74 


35 
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a Weighted average of 2.80 ± 0.38 g§ and 2.44 ± 0.26 |45| 
b Weighted average of 76 ± 8 @ and 72 ± 13 (47 1 
c Ref. @ 
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FIGURES 



FIG. 1. The Gamow peak for the Maxwellian and the Tsallis' distribution. The exponen- 
tially decreasing solid (dashed) curve labeled Maxwell (Tsallis) shows the Maxwellian (Tsallis with 
q = 0.98) energy distribution. Both the energy distributions have been normalized to a huge 400 
at E/kT = 0, to emphasize their tiny difference (the part of the distributions shown contains 
only about 0.3% of the total area). The rapidly increasing dotted curve shows the behavior of 
the penetration factor exp — yj Eg / E with Eq = 4000/cT corresponding to r = E$/kT = 10; the 
normalization is arbitrary. The solid (dashed) peak shows the product of Maxwellian (Tsallis) 
distribution function times the penetration factor. 

FIG. 2. The spectrum of solar neutrinos produced by the pp reaction. The solid line shows the 
spectrum at kT = 1.36 keV (about the solar central temperature) with the Maxwellian distribution 
for the proton relative energy, the dashed line the spectrum with the Tsallis' distribution at the same 
temperature. For comparison the dotted line show the spectrum with a Maxwellian distribution 
at T = 0. 



FIG. 3. The solar neutrino spectrum. For continuous sources, the differential flux is in 
cm _2 s _1 MeV -1 . For the lines, the total flux is in cm _2 s~ 1 . There are also indicated the thresholds 
of the three kinds of present experiments. The radiochemical experiment at Homestake [37,48] 
that uses using the reaction v e + 37 C1 — > 37 Ar (Chlorine) has the threshold 0.814 MeV. The 
radiochemical experiments Gallex 1 35 ,46] and SAGE |39|j47[] that use the reaction f e + 71 Ga — > 71 Ge 



(Gallium) have the threshold 0.2332 MeV. The thresholds of the Kamiokande and Superkamiokande 
experiments p5| , p4| , ^] , which detect the Cherencov light emitted by the scattered electrons in water 
(Water Cherenkov) are chosen between 5.5 and 8 MeV. 
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